A0-A113  499  NAVAL  RESEARCH  LAB  WASHINGTON  DC  F/6  20/9 

DRIFT  RESISTIVE  INTERCHANGE  AND  TEARING  MODES  IN  CYLINDRICAL  GE--ETC(U) 
APR  82  J  M  FINN»  W  M  MANHEIMER»  T  M  ANTONSEN 
UNCLASSIFIED  NRL-MR-4783  NL 


3499 


security  Classification  of  this  page  fWh*n  Dmtm  Emorod) 


REPORT  DOCUMENTATION  PAGE 


REPORT  n UMBER 


NRL  Memorandum  Report  4783 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.I  >■  RECIPIENT'S  CATALOG  NUMBER 

A $-A/f3 


«.  TITLE  (Mtd  Submit)  S  TYPE  OF  REPORT  A  PERIOD  COVEREO 

DRIFT  RESISTIVE  INTERCHANGE  AND  TEARING  MODES  Interim  report  on  a  continuing 
IN  CYLINDRICAL  GEOMETRY  NRL  problem. _ 

«.  PERFORMING  ORG  REPORT  NUMBER 


7.  AUTNORf*; 

J.M.  Finn,  W.M.  Manheimer,  and  T.M.  Antonsen,  Jr.41 


S.  CONTRACT  OR  GRANT  NUMBERO; 


f.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 


Naval  Research  Laboratory 
Washington,  DC  20375 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


Department  of  Energy 
Washington,  DC  20545 


12.  REPORT  DATE 

April  13, 1982 


is  NUMBER  OF  PAGES 

51 


4.  MONITORING  AGENCY  NAME  A  AODRESS (H  differ**!  from  Controlling  Oftlco)  »S-  SECURITY  CLASS,  (of  thio  ropofl) 

UNCLASSIFIED 


i«.  distribution  statement  t*i  ink  r.pomj 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  tho  obotroct  ontorod  In  Block  20,  ft  dllloronl  from  Rmport) 


l|.  supplementary  notes 

♦Present  address:  Science  Applications,  Inc.,  McLean,  VA  22102 
This  work  was  supported  by  the  Department  of  Energy. 

1 19.  KEY  WOROS  (Centlnw  on  r«v«r««  tIR  II  n«c«««A ry  Pid  Idwnttty  by  tlock  "numbme) 


Reversed  field  pinch 
Resistive  interchange  instability 
.Hall  terms  in  plasma  instability 


ecte; 

APR  1  5  OT 


tfl-^y^BSTRACT  (Conttnum  on  rororoo  of  do  ff  nocoooory  md  Idontlty  by  block  numfbot) 

Resistive  interchange  and  tearing  modes  are  studied  including  the  effects  of  electron  parallel  pressure 
gradient  and  ion  polarization  drift  (Hall  terms)  as  well  as  finite  plasma  compressibility  and  perpendicular 
resistivity.  It  is  found  that  for  unfavorable  curvature  and  weak  tearing  forces  (positive  but  small  &')  ail 
modes  are  stabilized  provided  the  effects  of  the  Hall  (drift)  terms  are  large  enough.  The  principal 
stabilizing  influences  are  plasma  compression  and  coupling  to  outgoing  drift  waves.  The  relevance  to 
reversed  field  pinch  and  spheromak  experiments  is  discussed. 


DD  |  JAN*!  M73  COITION  OF  I  NOV  «t  IS  OBSOLCTC 


S/N  0102-014- 4401 


sc curity  classification  of  this  raoc  (W*m  sm  Kmtw; 


'  'O  i 


CONTENTS 


I.  INTRODUCTION  .  1 

II.  THE  EQUATIONS  AND  SCALINGS  .  5 

III.  APPROXIMATE  ANALYTIC  SOLUTIONS  .  12 

IV.  NUMERICAL  RESULTS  .  25 

V.  CONCLUSIONS  .  39 

ACKNOWLEDGMENTS .  40 

REFERENCES  .  41 

APPENDIX .  43 


iii 


v-: 


f  masowi  p-os  ilanww  nww 


DRIFT  RESISTIVE  INTERCHANGE  AND  TEARING  MODES  IN  CYLINDRICAL  GEOMETRY 


I.  Introduction 

In  a  plasma  with  unfavorable  average  curvature  along  each  field  line,  it  has 
1-3 

long  been  appreciated  that  finite  conductivity  weakens  the  shear  stabilization 
criterion  of  ideal  interchange  modes  (the  Suydam  criterion^1  or  its  toroidal 
generalization,  the  Mercier  criterion).^’6  In  a  sheared  slab  with  density 
gradient  and  with  gravity  simulating  curvature,  and  with  incompressible 
displacements,  the  plasma  is  unstable  to  an  infinite  number  of  such  resistive 
interchange  (or  resistive  g)  modes  localized  near  each  mode  rational  surface 
(at  which  k  •  =  0,  where  k  is  the  wave  vector). 

In  a  cylindrical  model  with  curvature  and  pressure  gradient  these  results 

•i 

are  strongly  modified  by  plasma  compression.  In  a  plasma  with  very  low  pressure 
(but  finite  pressure  gradient)  the  results  are  formally  identical  to  those  in  the 
slab  model.  Finite  plasma  pressure  introduces  the  stabilizing  effect  of  plasma 
compression  which,  If  large  enough,  stabilizes  all  but  one  mode  for  each  mode 
rational  surface.  This  mode,  whose  magnetic  perturbation  is  not  totally 
localized  to  the  mode  rational  surface,  and  which  couples  to  the  exterior 
(magnetohydrodynamic)  region,  is  in  fact  the  tearing  mode.  For  unfavorable 
curvature  and  stabilizing  tearing  forces  from  the  outer  region  (A  <  0,  defined 
in  Sec.  Ill),  it  is  driven  unstable  by  curvature  in  the  same  manner  as  the 
localized  resistive  interchanges.  It  has  recently  been  pointed  out^  that  for 

a 

large  |A  |  or  small  D,  the  tearing-like  coupling  to  the  exterior  region  has  a 

strong  stabilizing  influence  that  reduces  the  growth  rate  from  the 
2/3 

conventional  D  scaling  (D  is  the  Suydam  parameter,  defined  in  Sec.  Ill)  to 
4  '4 

a  D  /A  scaling.  This  transition  has  also  been  observed  in  numerical 

D 

simulations . 

It  has  been  suggested  that  resistive  interchange  modes  are  responsible  for 

o 

transport  in  reversed  field  pinches,  which  have  unfavorable  average  curvature 
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because  their  safety  factor  q  (defined  in  Sec.  II)  is  less  than  unity.  It  is 

also  possible  that  fluctuations^  which  are  seen  late  in  the  discharge  in  the  Los 

Alamos  experiment  (ZT-40)  and  which  are  predominately  m  ■  0  (m  is  the  poloidal 

2/3 

mode  number),  are  resistive  interchange  modes.  The  transition  from  D  behavior 
4  '4 

to  D  /A  behavior  in  reversed  field  pinch  geometry  occurs  at  g  —  7  %  for 

3  6 

magnetic  Reynolds  number  R  *  10  but  at  0  ~  1  %  for  R  -  10  ,  so  that  present 

m  m 

day  devices  (with  g  ~  10%)  should  be  in  the  larger  growth  rate  (i.e. 

2/3 

D  )  regime.  We  shall  discuss  this  point  further  in  Sec.  IV. 

In  addition,  compact  torus  devices  such  as  the  spheromak,  which  have  q  <  1, 

and  which  therefore  have  unfavorable  average  curvature,  should  be  subject  to  such 

instabilities.  It  has  been  recently  found  that  spheromaks  are  optimally  stable 

to  ideal  magnetohydrodynamic  tilt  and  shift  modes  if  their  flux  surfaces  near  the 

11  12 

magretic  axis  are  nearly  circular.  *  Also,  transport  studies  indicate  that 
the  toroidal  current  and  temperature  become  quite  peaked  near  the  magnetic  axis, 
as  in  a  tokamak.  Therefore,  in  spite  of  the  small  toroidal  aspect  ratio  of 
spheromaks,  which  may  cause  resistive  instabilities  to  have  some  ballooning  like 
structure  (i.e.  to  be  sensitive  to  local  as  well  as  average  curvature),  the 
variation  of  the  field  line  curvature  along  field  lines  in  the  region  of  large 
pressure  gradient  may  be  small. 

The  purpose  of  the  present  paper  is  to  investigate  the  drift,  or  Hall 
effects,  on  these  instabilities.  These  terms  become  important  in  the  component 
of  Ohm's  law  parallel  to  the  magnetic  field  when  |u|  <  |w*el  (u>*e  is  the  electron 
diamagnetic  frequency  kcT^/eBL^,  and  L^  is  the  density  gradient  scale  length)  and 
represents  the  fact  that  a  parallel  component  of  the  electric  field  may  be 
balanced  by  parallel  electron  pressure  gradient  as  well  as  collisions.  The 
inclusion  of  these  terms  in  the  perpendicular  component  of  Ohm's  law  is 
equivalent  to  the  inclusion  of  the  ion  polarization  drift.  As  we  shall  see,  the 
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effect  of  the  Inclusion  of  the  Hall  terms  on  these  modes  is  to  impart  a  real  part 
to  their  frequency  <o  <  w.  and,  if  these  effects  are  strong  enough,  to  stabilize 

the  modes  altogether.  The  m  ■  0  fluctuations  in  ZT-40  have  a  frequency  which  is 

,  ^  10 
of  order  o>.  . 

*e 

The  model  we  use  has  cylindrical  geometry,  and  thus  has  field  line  curvature 
rather  than  a  fictitious  gravity.  The  only  toroidal  effect  is  that  the  axial 
wave  number  is  quantized  k  «  *  -  n/R,  where  R  is  the  major  radius.  We  include 

finite  plasma  compressibility  (with  electrons  isothermal  along  field  lines  and 
adiabatic  ions)  and  perpendicular  resistivity,  which  is  neglected  in  tokamak 
ordering. 


In  Sec.  II  we  explain  the  physical  model  in  more  detail  and  derive  the  basic 
equations. 

In  Sec.  Ill  we  derive  analytic  results  in  several  limiting  cases,  including 
the  incompressible  "slab"  limit  (p  ♦  0,  dp/dr  *  0)  where  plasma  compression  and 
perpendicular  resistivity  are  ignored.  We  recover  the  dispersion  relation  for 
the  drift  tearing  mode1**'15  from  a  more  general  dispersion  relation  for  drift¬ 
tearing-resistive  interchange  modes.  This  more  general  result  shows  that  slab¬ 
like  resistive  interchange  modes  remain  unstable  for  arbitrarily  large  /Or- 

where  Qf  is  a  nominal  resistive  instability  growth  rate  (defined  in  Sec.  II)  but 

-1/2 

with  a  growth  rate  that  scales  as  «  .  We  also  derive  a  dispersion  relation 

with  perpendicular  resistivity  and  plasma  compression  for  localized  modes 

with  mjke  large;  this  result  shows  that  for  large  no  such  mode  is  unstable. 

In  the  Appendix  we  generalize  this  to  less  localized  modes  which  couple  to  the 

outer  region  through  outward  propagating  drift  waves  and  show  that  no  instability 

+ 

exists  for  large  w*  /Q  unless  A  (defined  in  Sec.  Ill)  is  above  a  critical  value 
*e  r 

which  is  positive.  It  is  possible  that  the  fluctuations  seen  in  reversed  field 
pinches  are  resistive  Interchange  modes  which  become  unstable  when  the  mode 


rational  surface  is  near  the  edge  of  the  plasma,  where  compressional 
stabilization  is  at  a  minimum. 

In  Sec.  IV  we  show  numerical  results  which  are  consistent  with  the  analytic 
results  of  Sec.  Ill  in  showing  that  no  instability  exists  for  large  enough 
w*e/Qr,  when  the  coupling  to  outward  propagating  drift  waves  is  strong.  We  find 
the  critical  values  for  for  the  localized  and  nonlocalized  (tearing)  mode  and 
find  these  values  below  but  very  near  typical  values  for  reversed  field  pinch 
parameters. 


II.  The  Equations  and  Scalings 

In  this  section  we  review  the  derivation  of  the  equations  for  resistive 
interchange  modes  including  the  additional  effect  of  the  Hall  term  in  Ohm's 
law.  As  we  will  see,  these  terms  allow  for  the  coupling  of  resistive  interchange 
modes  to  drift  modes.  Since  it  is  now  known  that  drift  waves  are  stable  in  slab 
geometry  in  a  sheared  field,  this  coupling  is  a  stabilizing  effect. 

We  adopt  the  scaling  as  assumed  in  Coppi,  Greene  and  Johnson^  and  our 

previous  work.^  Outside  a  small  distance  1^.  away  from  the  rational  surface,  the 

ideal  MHD  equations  hold,  but  within  this  distance  Lf,  finite  resistivity, 

perpendicular  inertia,  and  Hall  effects  (electron  pressure  gradient  and 

polarization  drift)  are  important.  We  will  assume  that  L^/a  -  c,  where  a  is  the 

radius  of  the  mode  rational  surface  and  a  is  also  assumed  to  be  roughly  equal  to 

the  scale  of  variation  of  the  eigenfunction  in  the  other  direction  perpendicular 

to  B^.  Since  the  plasma  is  assumed  to  be  stable  to  localized  ideal 

magnetohydrodynamic  modes  (Suydam  stable),  the  growth  rate  also  scales  as  e.  In 

1/3 

this,  the  so-called  slow  interchange  ordering  of  Ref.  3,  e  also  scales  as  n  . 

In  cylindrical  geometry,  the  eigenfunction  has  the  form  f(r)  exp 
(yt  +  im0  +  ikz)  where,  in  a  large  aspect  ratio  torus,  k  ■  -n/R.  A 
superscript  ~  denotes  a  perturbed  quantity.  The  mode  rational  surface  is 
determined  by 

m  -  nq,  (1) 

where  q  *  rB  /RB  .  The  operator  B  •  7  operating  on  a  scalar  is  given  by 

Z  w  “  “ 

s 

inBeq  2 

II  •  V  - - - -  (r  -  a)  +  0(r  -  a)  (2) 

where  r  •  a  is  the  location  of  the  mode  rational  surface;  within  the  singular 
region  the  right  hand  side  of  (2)  is  of  order  e.  (We  use  this  form  rather  than 

5 


that  of  Ref •  3  since  it  applies  also  to  m  *  0  modes,  which  have  their  rational 

surface  where  B  *  0.) 

z 

Since  the  resistive  interchange  mode  is  characterized  by  nonzero  resistivity 
allowing  the  fluid  to  slip  through  the  sheared  field  near  the  rational  surface, 
the  mode  is  mostly  fluid  flow  with  only  a  small  magnetic  perturbation. 

Thus  jl  «■  e  where  _£  is  the  linearized  fluid  displacement.  Since  the  growth 
rate  is  assumed  to  be  slow  compared  to  the  magnetosonic  speed,  we  have 

B  •  B  +  p  -  0.  (3) 

where  p  is  the  perturbed  total  (electron  plus  ion)  pressure.  However,  for  Ohm's 
law,  we  use  the  full  ion  momentum  equation 

M  3v 

E  +  vxB-nJ=  7  "at  +  M- 

where  p^  is  the  ion  pressure,  p  is  the  density  and  M  is  the  ion  mass.  We  assume 
that  the  ions  obey  an  adiabatic  equation  of  state,  i.e.  the  perturbed  ion 
pressure,  pA  is  given  by 

Pi  -  -  pi  5r  ~  r  pi CV  •  £),  (5) 

where  T  is  the  specific  heat  ratio  and  £  is  the  fluid  displacement.  Since  both 
terms  on  the  right  hand  side  of  Eq.  (5)  contribute  roughly  equally  to  the 
perturbed  ion  pressure,  7  •  £  is  of  order  e. 

The  electrons,  on  the  other  hand,  are  assumed  to  be  isothermal  along  a  field 
line,  or 

B  •  VT+B  •  7T  -0, 

—  -  e  —  —  e 

6 

’  H 


which  gives 


a  B 

T  -  - -  T  ,  (6) 

e  e 

inB0q  x 

where  henceforth  x  =  r  -  a  and  the  perturbed  electron  pressure  is  pg  “ 

(pT  +  pTg)/M.  Note  that  if  Bf  has  even  symmetry  about  the  field  line,  is 

singular.  This  singularity  can  be  removed  by  accounting  for  finite  electron 

thermal  conduction  along  a  field  line.  However  if  has  odd  symmetry,  or 

if  T  ■  0  at  r  ■  a,  there  is  no  singularity.  If  we  define  a  local  coordinate 
e 

system  i  i  and  i,  *  i  x  i.  where  i  and  i,  are  unit  vectors  in  the  direction 

— r,  —  I  — 1  -r  —  I,  -r  —I 

of  r_  and  j^,  then 

I  *  Cr  ir  +  5I  -«  +  -i» 


B  1 
r  — i 


+  BJ  -I  +  B1  -1’ 


where  £  ~  e°,  C  B  ~  e  and  B  ~  e2.  Then  the  condition  that  7  •  £  ■  0  to 

II  r |  X) 1  r 

lowest  order,  J?  •  S  *  0,  pressure  balance  and  the  continuity  equation  gives  the 
following  four  equations 


1RBe  Hr 

+  nB  dr 


0, 


*  .  1RBe  3Br  . 

B1  +  nB  dr  “  °’ 


BB|  -  Pi  5r  +  rPl(V  •  £)  -  pT^/M  - 


aB 


re 


InB^q  xM 


Te  * 


P  ■  -  5fP  -  P(V  •  .£)• 


(7) 

(8) 

C9) 

(10) 
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Two  other  equations,  which  are  the  same  as  in  Ref.  3,  come  from  operating  on  the 

-2 

total  momentum  equation  with  B  •  and  with  either  7  •  (B  B  x)  or  B  •  7  x  are 


9  iftBfiq  Bx  ^ 

PY  «  -  <Pe  +  Pl)  Br - —  -  B,. 


2d\  2Z2*-  d\ 

PY  - T  -  — —  B.- 


dr 


R2a 


dr 


(ID 

(12) 


Equation  (12)  is  called  the  annihilated  momentum  equation,  since  the  operator 
annihilates  information  about  propagation  along  the  lines  of  force  and  also 
information  about  7(1$  •  B_  +  p).  In  order  to  neglect  collisional  and 
collisionless  ion  viscosity  in  Eq.  (12),  we  require  T  «  Tfi.  The  next  two 
equations  are  the  radial  and  parallel  components  of  Ohm's  law.  It  is  these  two 
equations  that  contain  the  effects  of  the  additional  Hall  terms.  They  are 


2~ 
d  B 

£  H _ r 

r  y  2 

dr 


‘"V  *  M  in.  _ 

a  *r  e  '  R.B  V 


(13) 


2~ 

d  B  -inBaq  x  ,  A  D2  ' 

B  -  — - r- - e——  l  -  B(7  •  O  -  % 

I  Y  J  2  a  v-  B  ‘’r 

dr 


Ml  .  _M_  InB 
e  3r  +  p2e  B0Ry 


(P 


~  d£  . 


i  dr 


(14) 


As  in  Ref.  3  the  eight  equations,  (7)  -  (14),  can  be  reduced  to  three  second 
order  differential  equations  for  the  three  variables  B^,  B (  and  Cr*  We  introduce 
the  variables  used  in  Ref.  3:  Q  ■  y/Qj.*  C  ■  Cf »  X  ■  x/Lf,  and 


¥ 


-  flr » 
”“8"  Lr 


T 


-2Ba 


R2q 


2B2 

e 


9 


8 


n 


\ 


2p  a 


2  2  "2  ’ 

Be  R  * 


B  =  2p(a)/B  , 


2  '2 

S  =  4/R  q  % 


(2  2  \  1/6 

-VILA _ I 

2_22l 

0  q  / 


Or 


(2_  2  '2  \  1/3 

^  Be  q  \ 

? )  ' 


as  well  as  two  new  variables  g  (to  measure  the  strength  of  the  Hall  effect) 


R  ,/3,»,'.,‘/3  V 

8  -  J - = - <^T>- 

I Rq  I  ci 


(15) 


where  <m  .  is  the  ion  cyclotron  frequency  using  B  ,  R  is  the  magnetic  Reynolds 
ci  „  o  m 

aV  - 

number,  R  =  — -  ,  V.  is  the  Alfven  speed  using  B  ,  and  k  *  (1  -  Tp  p  / 

m  h  A  0  1 

pi  p^pi  ^p  In  terms  these  variables  the  three  coupled  second  order 
equations  reduce  to 


d2S  j_  »  dff  1  T 

dX2  Q2  dX2  Q2 


Q(f  +  X  5)  +  IQ*  »  +  X  , 


(16) 

(17) 
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(IB) 

—  -Q(l+  \  +  pV  T+Q  (S-D  -  ~)  C+  ^  XT  +  IQ*  ~  <(T-DO-2iQ2g 

dX  Q  cp  c  ^  c  dX 

where  Qt  *  2gD/S. 

In  deriving  Eqs.  (16)  -  (18),  we  have  further  assumed  that  7Tg  «  0  so  that 

subtleties  concerning  singular  behavior  of  electron  temperature  are  avoided. 

The  T  in  Eq.  (18)  is  then  a  combined  specific  heat  ratio  which  accounts  for  the 
c 

fact  that  the  electrons  are  isothermal  but  the  ions  are  adiabatic. 


r 

c 


pTe/M  +  rPi 

p  +  p, 
re  i 


(19) 


If  the  ion  temperature  is  zero,  Tc  *  1. 

The  terms  in  Eqs.  (17)  and  (18)  which  contain  the  g  factor  are  the  effects 
of  the  Hall  terms.  Since  they  are  all  purely  imaginary,  they  provide  a  real  part 
of  the  frequency.  Notice  that  the  terms  in  Eq.  (18)  proportional  to  ic  are  zero 
for  a  "baroclynic"  ion  pressure  profile,  that  is 


1  d£ _ 1_ 

P  dx  '  Tpi 


(20) 


r 

or  Pj  '  p  .  In  our  numerical  computations,  we  generally  make  this  assumption. 
Also,  we  have  neglected  collisionless  and  collisional  ion  viscosity  and  thermal 
conduction.  For  reversed  field  pinch  parameters,  this  can  only  be  justified 
if  Tg  »  Tj.  Finally,  to  make  contact  with  conventional  drift  waves  and  drift 

A 

ordering,  let  us  note  that  the  drift  wave  frequency,  mcTgp  /eBap  is  given  by 


u 


*e 


Q*V 


(21) 
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As  we  will  see  in  the  next  section,  'combinations  like  Q  +  IQ*  (proportional 
to  a)  -  m*e)  frequently  appear  in  the  analysis  of  Eqs.  (16)  -  (18)  and  it  is 
through  these  terms  that  resistive  interchange  and  tearing  modes  couple  to  stable 

drift  waves. 
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III.  Approximate  Analytic  Solutions 

There  are  at  least  three  limits  in  which  one  can  get  analytic  solutions  to 

Eqs.  (16)  -  (18).  The  first  is  g  =  0  and  these  have  been  discussed  fully  in 

Refs.  1  -  3,  7.  The  second  is  the  unphysical,  yet  illuminating,  limit  Tc  +  0  and 

the  third  is  the  limit  T  =  0,  g30  »  1. 

We  begin  by  discussing  the  limit  T ♦  0.  In  this  limit,  plasma  compression 

makes  no  contribution  to  perturbed  pressure.  In  this  sense,  it  is  rather  like 

the  incompressible  theory  of  resistive  interchange  modes  in  a  slab,  where  the 

curvature  -  pressure  gradient  term  in  Eq.  (12)  (proportional  to  )  is  replaced 

by  a  gravity  -  perturbed  density  term.  As  we  shall  see,  the  results  of  this 

limiting  case  are  similar  to  the  simple  slab  mode  without  Hall  terms,  in  the 

sense  that  an  unstable  mode  exists  for  arbitrarily  large  shear  and  g. 

If  T  =0,  then  Eq.  (18)  shows  that  T  *  D  £,  so  that  Eqs.  (17)  and  (16)  become, 
c 

respectively 


~  -  (Q  +  iQ*)  (y+X()  (22) 

dX 

2  2 

4-4  -  K  (Q  +  IQ*)  ?  -  K  +  K  (Q  +  iQ*>  f  (23) 

dx^  <r  Q  Q 

2  2 

where  we  have  used  Eq.  (22)  to  eliminate  the  d  T/dX  term  in  Eq.  (16).  To  solve 

Eqs.  (22)  and  (23)  we  first  make  the  constant  Y  approximation.  If  £  is  even  in 

X,  ¥  is  odd,  so  that  with  this  approximation  f  -  0.  In  this  case  the  lowest 

2 

order  eigenfunction  £  is  given  by  £  ■  exp(-  o  X  ),  where 


o 


i(S 

2  v 


+  iQ* 


2 


(24) 


and  the  eigenvalue  (i.e.  dispersion  relation)  by 
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12 


(25) 


1 

Q(Q  +  i  Q*)2  -  D. 

The  eigenfunction  is  only  valid  if  it  approaches  zero  as  X  +  00 .  Using  Eq. 
(25)  and  making  use  of  the  fact  that  D  is  positive  real,  we  find  that  the 
eigenfunction  is  well  behaved  for  X  ♦  +  «  [Re(a)  >  0]  as  long  as 

Re(Q)2  >  Im(Q)2. 

If  g  =  0,  Eq.  (25)  is  the  standard  dispersion  relation  for  resistive 
interchange  modes.  As  g  increases,  the  real  part  of  Q  decreases  while  the 
Imaginary  part  increases.  In  the  limit  of  g  + 

_ 

Q  -  D(iQ*)  2 

so  that  the  condition  for  validity  of  the  eigenfunction  becomes  marginally 
satisfied  for  g  <»;  note  also  that  Q  0  in  this  limit.  Thus  if  **  0,  the 
theory  predicts  that  the  Hall  effect  reduces  the  growth  rate  of  a  resistive 
interchange  mode,  but  does  not  stabilize  the  system  for  finite  g.  These  results 
are  summarized  in  Fig.  1. 

It  is  Just  as  easily  shown  that  the  higher  order  eigenfunctions  have 

2 

eigenfunction  C  =  H^(z)e  Z  ^  for  even  4,  where  z  *  [(Q  +  iQ#)^/Q^2]  X,  and 
are  Hermite  polynomials.  The  eigenvalues  for  these  modes  satisfy 
Q(Q  +  1Q*)1/2  -  D/<24  +  1). 

The  remaining  thing  to  check  is  the  validity  of  the  constant  ¥ 
approximation.  As  we  will  see,  it  is  valid  if  D  «  1,  which  holds  for  a 
low  3  plasma.  According  to  Eq.  (24),  the  length  scale  of  the  eigenmode 
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2  ' 

is  [(Q  +  IQ*)  /Q  ]  •  Using  this,  we  find  that  the  Y  terms  on  the  left  and 

1/2 

right  hand  side  of  Eq.  (22)  are  related  as  1  :  (Q  +  iQ^)  Q.  Thus  the  left 

2  2 

hand  side  of  (23)  dominates  for  D  «  1,  ^  ving  (to  lowest  order)  d  Y/dX  *  0. 

Thus  if  r  *  0,  it  is  possible  to  derive  the  dispersion  relation  and 

eigenfunction  for  all  resistive  interchange  modes  in  a  cylindrical  plasma.  For 

the  case  of  g  *  0,  they  are  identical  to  resistive  interchange  modes  in  a  slab, 

when  an  artificial  gravity  is  used  in  place  of  curvature.  As  g  increases,  the 

frequency  behaves  as  in  Fig.  1  (for  all  even  l  modes.)  Stability  is  not  achieved 

-  1/2 

for  any  finite  g,  but  the  growth  rate  behaves  as  g  as  g  ♦  with 

|lm/Qi/|ReQ|  approaching  unity  as  g  The  latter  observation  implies 

that  Re (a)  approaches  zero  as  g  +  «  so  that  the  outgoing  drift  waves  are  not 
spatially  damped  for  g  +  ®. 

We  now  consider  modes  with  odd  £  and  even  Y  which  couple  to  the  outer  MHO 

region.  The  problem  now  is  to  solve  Eq.  (23)  in  which  Y  is  regarded  as  a 

constant.  Since  we  are  looking  for  modes  which  do  not  satisfy  the  homogeneous 

equation  for  £  (these  are  the  resistive  interchange  modes  we  have  just 

discussed),  we  are  interested  in  the  particular  solution.  This  can  be  found  most 

ikX 

easily  by  Fourier  transforming  the  equation.  Multiplying  by  e  and  integrating 
over  X,  we  find 


Q  +  iQ- 


dk 


2"  D 

+  k^£  -  £ 

Q 


Q  +  iQ* 

2iri  - y ~  *0  6  <k> 


(26) 


where  we  have  used  the  notation  £(k)  ■  /  "^dX  £(X)e^^.  Thus  we  want  to  solve 

A 

Eq.  (26)  subject  to  the  boundary  condition  that  £  is  well  behaved  as  k  *  +  «. 
Since  the  inhomogeneity  is  the  derivative  of  a  delta  function,  we  simply  solve 
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the  homogeneous  equation  for  |X|  >  0  and  calculate  the  discontinuity 
across  X  *=  0.  For  k  >  0,  Eq.  (26)  is  the  equation  for  parabolic  cylinder 
function.  For  positive  k  we  find 


where 


(27) 


(28) 


(29) 


and  A  is  a  coefficient  to  be  determined  shortly.  The  function  U(a,ic)  is  the 
parabolic  cylinder  function  (in  the  notation  of  Ref.  17)  which  approaches  zero 

as  k  approaches  infinity  along  the  real  axis.  For  a  conventional  tearing 

* 

mode,  Q*  «  0  and  Q  is  real,  so  £  does  have  the  proper  behavior.  If  the  Q*  terms 
are  included,  the  situation,  as  we  will  see,  is  somewhat  more  complex.  We  oust 
check,  a-posteriori  ,  that  the  eigenfuctlon  derived  by  the  dispersion 
relation,  U(a,<)  has  the  proper  behavior  as  k  +  ». 

A 

Equation  (27)  is  the  functional  form  of  £(k)  for  positive  k.  Since  the 

A 

source  term  on  the  right  hand  side  of  Eq.  (26)  has  odd  symmetry  in  k,  £(k)  oust 
also  be  an  odd  function  for  the  particular  solution.  Thus  we  find 


£(k)  -  -  A  U(a,-  k) 


k  <  0 


(30) 


The  value  of  A  can  be  determined  by  integrating  Eq.  (26)  across  the  derivative  of 
the  delta  function  at  k  *  0.  The  result  is 
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A 


-  Hi  » 

_ £ 

U(a,o) 


(31) 


To  summarize,  the  solution  to  the  inhomogeneous  equation  for  £00 >  Eq.  (26),  is 
given  by  Eqs.  (27)  -  (31). 

To  get  the  dispersion  relation,  we  assume,  as  usual,  that  the  outer  ideal 

* 

magnetohydrodynamic  region  specifies  a  value  of  A  ,  the  discontinuity 
in  (dB^/dO/B^,  across  the  singular  layer.  Matching  with  its  counterpart  from 
the  inner  region,  we  obtain,  from  Eq.  (23) 

L.A*  -  JIdX(Q  +  iQ*)(Y  +  X£)  (32) 

1  rm  A  “iVcX 

To  do  the  integrals  over  X  note  that  £(31)  «  —  J_co£(k.)e  dk  while  S'  (assumed 
00  — ^  ^ 

constant)  is  ¥  j  a)6(k)e  dk.  Note  that  in  doing  the  £  integral  over  X,  one 
will  ultimately  have  a  derivative  of  a  delta  function  in  the  k  integral. 

Inserting  these  forms  into  Eq.  (32)  we  see  that  both  the  V  and  £  contributions 
have  a  delta  function  component  at  k  *  0,  but  these  two  divergent  contributions 
cancel.  This  cancellation  is  of  course  the  analog  of  the  separate  Integrals  in  X 
space  of  f  and  X  £  diverging  for  X  ♦  «,  while  the  integral  of  ¥  +  X  £  converges. 

The  remaining  part  of  the  k  integral  comes  from  the  derivative  of  the 
parabolic  cylinder  function  at  k  ■  0.  The  final  result  is 

(  —  —  )  1  i 

L  &  -  2  *  F  -  4-+  l -  (Q  +  iQ*)4  Q2  (33) 

r<±  +  f) 

If  0  and  Qa  are  both  equal  to  zero,  this  is  the  standard  dispersion  relation  for 
the  tearing  mode.  To  get  the  growth  rate  with  *  0,  the  procedure  is  to  solve 
Eq.  (33)  and  remain  on  the  branch  which  gives  the  tearing  mode  as  ♦  0.  For 
D  -  a  ■  0,  this  gives  the  standard  drift  tearing  dispersion  relation*4' ^ 
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CQ  +  *  w*iere  ^  is  ttie  Pure  (Q*  m  0)  tearing  mode  growth 

rate.  However  we  must  also  check  that  this  branch  is  well  behaved  as 
X  or  k  -*■  ®.  This  requires  that  the  real  part  of  a  [c.f.  Eq.  (24)]  remains 
positive.  Inserting  for  Q  +  iQA  from  the  dispersion  relation,  with  D  =  0,  we 
find  that  this  means 


♦  > 


3tt 

8 


(34) 


where  tan(4>)  *  Re(Q)/im(Q).  As  shown  in  Fig.  2,  the  dispersion  relation  does 

approach  this  point  while  the  mode  is  unstable1^,  so  that  at  this  point,  the 

instability  abruptly  disappears,  or  possibly  shifts  to  another  branch  of  the 

dispersion  relation.  However,  we  also  show  in  the  next  section  that  if 
2  2 

the  d  T/dX  terms  are  included,  the  growth  rate  can  be  brought  to  zero  by 

increasing  while  the  eigenfunction  remains  well  behaved. 

For  D  >  0,  the  dispersion  relation  to  Eq.  (33)  gives  other  solutions  of  a 

more  interchange  -  like  nature.  For  example,  for  A  ♦  «  we  obtain  roots 

1/2 

associated  with  poles  of  the  T  function  in  the  numerator,  Q(Q  +  iQ*)  * 

1/2 

D/ ( 21  +  1),  where  l  is  now  odd.  For  A'  +  0  we  obtain  Q(Q  +  iQ^)  - 
D/(2£  -  1),  again  only  for  odd  Z. 

We  now  continue  with  our  examination  of  analytic  solutions  of  Eqs.  (16)  - 

(18).  We  first  make  the  constant  ¥  approximation,  and  begin  by  looking  at 

homogeneous  solutions  to  Eqs.  (16)  and  (18)  [i.e.  with  ¥  -  0]  in  the  limit  of 

2 

large  X.  Keeping  only  second  derivatives  and  terms  proportional  to  X  we  find 
that  the  equations  reduce  to 


fdi  .  id  r  +  21&  x2 

dX*  v  S(T 
18 


(35) 


•  ^ 


(36) 


-  £  t  -  2iQ2e  e. 

dX  ^ 

2 

In  the  limit  of  large  X,  we  assume  that  spatial  variations  are  as  exp  (-  oX  ). 
Using  this  spatial  dependence  and  taking  X  °»,  Eqs.  (35)  and  (36)  reduce  to 

‘•°2  ■  T  +  7$  To  (37) 

bvj 

T 

4  o2T  «  -  8i  Q2go>2£  (38) 

o  H  o 

so  that  we  obtain  four  solutions 

2  — 

«  *  ±  I  2/S  ^2^  C  S“+Q  >2  1  (39) 

g2  1 

The  two  expressions  for  a  are  quite  different  in  magnitude  as  long  as  »  — . 

b  V 

In  the  limit  of  large  g,  we  expect  the  frequency  to  be  roughly  a  drift  wave 
frequency  |Q  Q*  ~  8g,  so  that  the  roots  are  well  separated  if  S  ~  1  and 


g3e  »  1 


(40) 


We  will  assume  Eq.  (40)  is  satisfied,  and  also  that  Q  ~  8g  and  S  -  1  and 
examine  analytic  solutions  of  Eqs.  (16)  -  (18)  which  do  not  couple  to  the  outer 
MHD  region,  l.e.  for  even  £,  T  and  odd  Y.  Taking  the  negative  sign  inside  the 
bracket  in  Eq.  (39),  we  find 


a 


(41) 


Notice  that  if  Q  is  purely  imaginary,  as  it  nearly  is  for  a  drift  wave,  to  lowest 

20 


order  the  exponential  does  not  decay  in  space,  but  oscillates  more  and  more 


rapidly. 


( 


However  on  keeping  next  order  terms  in  the  expansion  of 
1/2 

)  ,  one  can  easily  show  that  a  has  a  real  part  even  if  Q  is  purely 


imaginary.  Thus  it  is  possible  to  have  spatially  localized  shear  damped  drift 


waves  (i.e.  real  part  of  o  >  0  and  R  (Q)  <  0. )  In  the  next  section  we  show  that 

e 

this  is  indeed  the  case. 


To  continue,  we  will  assume  lengths  are  scaled  by  the  o  given  in  Eq.  (41) 

1/2  1/2  ci  ""1  -1/2 

with  Q  ~  gf}  so  that  X  ~  (gQ)  '  ~  gB  '  and  g  B  .  Then,  comparing  the 

3 

terms  in  Eqs.  (16)  and  (17)  and  using  1/g  B  as  a  small  parameter,  we  find  that 
(16)  and  (17)  reduce  to 


s  x 


(42) 


a  simple  algebraic  relation  between  £  and  T.  Using  this,  and  using  the  same 
ordering  on  Eq.  (18)  we  find  that  it  reduces  to  a  Weber  equation 


2  2 

4Q«  d  T 


+  1  fs  +  (s  -  If)  (  »  T 


2 

where  we  have  asumed  T  «  0,  and  also  assumed  —  »  1.  If  S  were  equal  to  zero, 

I  ip 

the  term  in  the  square  brackets  would  be  proportional  to  the  drift  wave 
dispersion  relation  Q  +  iQ^«  0  where  •  The  effect  of  the  S  term  is  to 

lower  the  drift  wave  frequency  for  D  >  0  (that  is  for  "bad"  curvature).  With 
the  d^T/dX^  and  X2T  terms  included,  Eq.  (43)  gives  the  dispersion  relation  for 
shear  stabilized  drift  waves.  Assuming  Q  ■  -iu>  =  -  1Q.+  iTflg  in  the 
d^T/dX^  and  X^T  terms  and  outgoing  wave  boundary  conditions,  we  find  the 
dispersion  relation  is 


21 


S  +  (S  -  |~)  C^)  -  -  2g(2£  +  1 ) / S 


re 


r  e 

c 


(44) 


-1  '2 

for  l  even.  If  the  shear  length  is  taken  as  =  |aq  /Rq  |  and 


l?vl 


l-1,  r 

n  c 


1  for  isothermal  electrons, 
2 


,  Mol  .  (2i  ♦  1) 

Q*  |u,*e<  L  B2 

s 


(45) 


For  a  slab  plasma  in  which  *  B,  the  damping  rate  is  the  classical  result  for 

1  R 

shear  damping  of  a  fluid  drift  wave  in  slab  geometry.  °  Thus  in  the  limit  of 
3 

g  B  »  1,  one  of  the  even  £,  T,  odd  ¥  modes  is  the  shear  stabilized  (damped) 
drift  wave. 

The  condition  for  the  validity  of  the  constant  ¥  approximation  for  this 

ordering,  i.e.  the  condition  for  the  dominance  of  the  left  hand  side  of  Eq.  (17), 

3  2 

is  easily  seen  to  be  g  0  «  1. 

We  will  continue  by  considering  the  other  choice  for  spatial  scale  length  in 
Eq.  (39), 


a 


.  8_ 

✓S 


(46) 


which  corresponds  to  a  much  more  spatially  localized  mode  than  the  shear 
stabilized  drift  wave.  Assuming  d/dX  ~  g^2  /s*^\  X  ~  S^Vg^2, 

3 

Q  -  0g  and  g  8  »  1,  we  once  again  examine  Eqs.  (16)  and  (18)  (assuming 
again  T  •  0;  In  this  ordering  the  constant  ¥  approximation  requires  8  «  1)  and 
pick  out  the  dominant  terms.  They  are 


A  .  21&X?  T  _  J_  T 

dX2  SQ2  Q2 


(47) 
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*L  T_ 


rc6 


21 Q 


2  d^C 

dX2 


(48) 


Combining  these  to  form  a  single  equation  for  T,  we  again  arrive  at  a  Weber 
equation 


4^.^T  +  (|^+21g)T.o,  <*9) 

d£  c p 

the  eigenvalues  of  which  are 

Q  =  -  iT3g  -  rc0|g|  (24  +  1)/S2  (50) 

for  even  4.  Again  the  result  is  a  damped  mode  with  a  real  frequency  of  order  the 
drift  frequency. 

Thus  as  the  plasma  becomes  more  and  more  collisionless  (that  is  as 
g  -*■  “)  all  localized  (i.e.  even  £,  T,  odd  Y)  unstable  modes  disappear.  Hence  a 
localized  resistive  interchange  mode  which  is  unstable  for  g  =  0  must  be  damped 
as  g  <»  if  m  remains  a  localized  mode.  This  is  an  interesting  and  useful 
result  which  is  not  so  for  the  simplest  case  of  Tc  =  0.  In  that  case,  as  we 
showed  earlier  in  this  section,  the  unstable  resistive  interchange  mode  persists 
as  g  +  «.  Thus  it  is  not  only  the  Hall  effect,  but  also  plasma  compressibility 
which  needed  to  stabilize  localized  resistive  interchange  modes  in  a  cylindrical 
plasma. 

In  the  next  section  we  show  by  a  numerical  study  of  Eqs.  (16)  -  (18)  that 
localized  resistive  Interchange  modes  are  Indeed  stabilized  as  g  ♦  ••  However 
this  conclusion  seems  to  be  even  more  general.  We  find  that  both  tearing  modes 
with  positive  A  and  resistive  Interchange  modes  which  couple  to  the  outer  ideal 
magnetohydrodynamic  region  (but  have  A  <  0)  are  all  stabilized  by  increasing  g. 
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Furthermore,  as  shown  in  the  Appendix 
for  positive  A  if  g  is  large  enough 


-U!*U  -'■JAJ'JSWfl 


IV.  Numerical  Results 

In  this  section  we  describe  results  obtained  by  integrating  Eqs.  (16)  -  (18) 

numerically.  We  integrate  the  full  sixth  order  set,  because  the  conditions  for 

the  validity  of  the  constant  -  V  approximation  cannot  be  assured  to  hold  for  all 

3  2 

values  of  g.  Indeed,  in  the  drift  wave  region  this  criterion,  g  8  «  1,  holds 

only  marginally  for  values  of  interest.  Nevertheless,  we  integrate  only  the 
singular  layer  equations,  using  a  higher  order  version  of  the  A  method  usually 
used  in  conjunction  with  the  constant  -  V  approximation.  Specifically,  we  obtain 
from  (16)  -  (18)  the  first  two  terms  in  the  asymptotic  series  relating  £  and  T 
to  4*  at  large  X: 


[1  +  (D/Q  +  2igQ)/X2]  V/X, 

(51) 

♦  -  [i  -  sq2/dx2]  DV/X. 

(52) 

Using  these  expressions  in  the  annihilated  equation  of  motion,  we  compute  the 
symmetric  and  antisymmetric  components  of  f,  which  for  large  positive  X  are  given 
by 


V  -*■  f  [  1  +  DlnX  +  D2*nX 
s  o 

+  D2(inX)2/2J  +  (S-2)  ¥  Q2/6X2, 

o 

V  ♦V  X  [  1  -  DlnX  -  D2JtnX 
a  o 

+  D2(£nX)2/2]  +  Sfo  Q2/2X.  (53) 

The  first  terms  in  and  V  are,  respectively,  the  expansion  for  small  D  of  the 
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exterior  Suydam  solutions  (for  D  <  j-)  f  X^l  and  V  X**2  where  q.  „  = 

4  o  o  1,2 

1/2  ' 

(1  T  (1-4D)  ]/2.  We  match  t  /f  to  the  outer  (ideal  magnetohydrodynamlc)  region 

/ 

by  means  of  two  parameters  6^  ^  =  b¥o  !&^Q  for  X  -*■  +  «•  and  6^  for  X  ♦  -  <»), 
where  4*  is  assumed  to  be  given  by  V  -  a'V  +  M1  .  The  advantage  of  taking  the 

S  £1 

higher  order  terms  in  the  asymptotic  expansion  in  (53)  is  that  we  can  then 

perform  the  numerical  integrations  over  a  relatively  small  interval  in  X. 

We  first  discuss  numerical  solutions  which  correspond  to  the  analytic 

solutions  of  Refs.  3,  7  (g  -  0)  and  Sec.  Ill  (  T  *  0). 

c 

Figure  3  shows  the  results  obtained  by  this  code  for  the  odd  £  tearing 

resistive  interchange  with  g  -  0  and  A  <  0.  For  very  small  D,  we  see  that 

'  4 

the  Q  =  (2.12irD/4L  A  )  scaling  of  Ref.  7  holds,  showing  a  much  smaller  growth 

-  2/3 

rate  than  the  conventional  (A  *  0)  value  Q  =  (itD/4)  .  For  larger  D  the 

assumption  2D/TB  «  1,  in  which  compressional  stabilization  is  a  dominant  effect, 
breaks  down  and  the  growth  rates  found  by  the  code  are  larger  than  those  found  in 
Ref.  7.  Nevertheless,  for  reasonable  values  of  the  Suydam  parameter  D,  the 

2/3 

growth  rate  is  still  considerably  smaller  than  the  conventional  value  (itD/4) 

By  increasing  D  while  keeping  all  other  parameters  fixed,  we  pass  through  five 

4 

distinct  regions,  the  (D/A")  regime  of  Ref.  7,  a  transition  regime  due  to  the 

2/3 

diminished  stabilizing  effect  of  plasma  compression,  the  classical  D  regime,  a 
regime  where  the  constant  Y  approximation  breaks  down  and  finally  D  1/4,  where 
the  mode  becomes  unstable  in  ideal  magnetohydrodynamics.  For  realistic 
parameters  (e.g.  for  a  reversed  field  pinch)  these  regimes  overlap  so  strongly 
that  the  individual  regimes  cannot  be  unequivocally  recognized  in  the  numerical 
results. 

* 

We  have  also  used  the  code  with  -  .02,  L.^6^  -  -  .02  (LrA  -  +  .04),  TB 

-5  -3 

-10  ,0-10  and  S  -  0.1  to  investigate  the  low  beta  drift  tearing  mode 

discussed  in  Sec.  111.  We  find  excellent  agreement  in  growth  rate  and  in  real 


Growth  rate  of  even  f,  odd  (,  T  resistive  interchange  mode  as  a  function 
of  the  Suydam  parameter  D,  with  LA  ■  -  >06,  S  -  0.1,  r  {  ■  0.1. 


part  of  frequency  up  to  the  point  e  =  w  /y  “  1.89  (y  is  the  g  =  0  tearing  mode 

*  t  t 

growth  rate)  or  g  ■  4.21,  where  there  is  no  spatial  damping  as  X  ♦  ®,  i.e.  a  [of 

eq.  (24)1  is  pure  imaginary.  These  results  are  shown  in  Fig.  2. 

We  have  also  studied  a  low  rB,  &  -  0,  even  £  resistive  interchange  and 

found  that  the  numerical  results  were  in  excellent  agreement  with  (25),  which 

shows  that  the  unstable  resistive  interchange  approaches  marginal  stability 

as  g  -*■  “  (  ♦  °>),  and  that  Re(a)  [of  Eq.  (24)1  approaches  zero,  i.e., 

6  =  Arg  y  approaches  -x/4,  as  g  +  •.  These  results  are  summarized  in  Fig.  1. 

We  now  turn  to  a  study  of  drift  resistive  interchange  and  drift  tearing 

modes  in  a  cylindrical  plasma  with  parameters  relevant  to  current  reversed  field 

pinch  experiments.  Nominal  values  of  density,  temperature,  and  magnetic  field 

of  n  «  2  x  lO^cm  T  -  T,  ■  100  eV,  and  B  *  4  kg  are  used.  We  fix  TB  ■ 

e  i 

0.167  in  the  code,  consistent  with  the  above  parameters.  We  have  observed,  as 

first  noted  by  Greene,1®  that  for  fairly  high  beta  (more  than  a  few  percent)  and 

fixed  D,  the  stability  results  are  relatively  insensitive  to  B  (and  therefore 

insensitive  to  the  difference  between  an  isothermal  and  an  adiabatic  equation  of 

state).  Assuming  Bessel  function  profiles  B.  ■  J.(pr),  B  ■»  J  (pr),  with 

0  1  z  o 

2.4  <  pr  <  3.1  (r  is  the  radius  of  the  conducting  wall),  we  find  q  R  - 
w  w 

-  ya  "  -  2.4,  is  again  the  radius  of  the  mode  rational  surface.  We  have  assumed 

that  a  is  at  the  field  null,  as  is  appropriate  for  m  *  0  modes.  We  obtain 

S  *  0.7.  For  stability  to  all  localized  resistive  interchanges  with  g  •  0,  we 

pick  S  -  2D/TB  ■  -  1,  or  D  ■  0.14.  As  discussed  in  Ref.  7,  this  condition  gives 

reasonable  pressure  profiles  with  a  finite  but  very  small  pressure  at  the  wall. 

A  typical  value  for  R  ,  the  magnetic  Reynolds'  number,  for  radial  scale  lengths 

m 

6  ** 

5-10  cm,  is  R  ■  10  ,  and  for  2.4  <  pr  <  3.1,  A  r  -  -5  for  m  ■  0 
m  w  w 
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and  m  *  1  modes  (the  latter  with  their  mode  rational  surface  inside  the  field 


null).  We  find  ~  -  .06,  and  g  ■  5. 

In  Fig.  4  we  show  the  complex  growth  rate  Q  as  a  function  of  g,  for  the 

odd  £,  T,  even  f  mode,  with  all  the  other  parameters  given  in  the  previous 

paragraph.  We  also  show  the  line  id  *  or  Q  *  -  iQ^.  We  see  that 

<j)/b)A  -  1/2  throughout  the  unstable  range  and  that  this  mode  is  stabilized 

for  g  >  2.35,  well  below  the  nominal  value  of  this  parameter  for  reversed  field 

pinches.  The  mode  structure  for  various  values  of  g  is  shown  in  Figs.  5-8. 

For  g  ■  0.5,  as  seen  in  Fig.  5,  £,  T,  and  ¥  have  real  parts  resembling  g  -  0 

resistive  interchanges,  but  also  have  imaginary  parts.  For  g  *  2.0,  2.25,  2.40 

(Figs.  6-8),  one  can  see  evidence  of  the  localized  behavior  in  T  due  to 

perpendicular  resistivity  and  also  evidence  of  outgoing  propagating  drift 

waves.  (It  is  easily  seen  in  Figs.  6-8  that  the  radial  phase  velocity  is  inward, 

1 8 

and  we  know  from  the  local  dispersion  for  drift  waves  in  a  sheared  geometry 
that  the  radial  group  velocity  is  in  the  opposite  direction,  corresponding  to  an 
outgoing  wave  boundary  condition.) 

We  have  also  studied  the  even  £  resistive  interchange  with  smallest  radial 

mode  number  l  »  0  (i.e.  no  nodes  for  *  0).  We  pick  all  parameters  the  same  as 

for  the  odd  £  mode,  except  TB  *  0.0512,  so  that  with  g  »  0  the  two  modes  have  the 

same  growth  rate.  (By  decreasing  TB  without  changing  D  we  decrease  the 

stabilizing  Influence  of  plasma  compression.)  This  mode  is  then  stabilized  for  g 

-  4.0,  as  shown  in  Fig.  9.  This  larger  value  is  required  presumably  because  this 

more  localized  mode  is  more  weakly  coupled  to  outgoing  drift  waves.  For  small  g 

we  observe  to  «  u./3.  The  mode  structure  at  marginal  stability  is  shown  in  Fig. 
r  * 

10.  Again  evidence  for  localized  behavior  due  to  perpendicular  resistivity  and 
drift  wave  behavior  at  larger  X  is  present,  and  the  apparent  radial  phase 
velocity  is  negative,  indicating  an  outward  propagation  of  energy. 
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We  have  also  Investigated  cases  in  which  k  *  0,  i.e.  when  the  ion 
temperature  gradient  is  smaller  (<  <  0)  or  larger  (ic  >  0)  than  that  of  the 

r  -1 

baroclynic  profile  ~  p  .We  find  that  increasing  ion  temperature  gradient 
is  a  stabilizing  influence.  However,  these  results  do  not  necessarily  have 
immediate  implications  for  devices  with  ~  because  we  have  neglected  the 
collisionless  ion  viscosity  and  thermal  conduction  terms,  which  are  of  the  same 
order  of  magnitude. 

We  have  also  investigated  the  odd  £  mode  of  Figs.  4-8  for  A  >0,  where  it 
goes  over  to  a  tearing  mode.  As  seen  in  Fig.  11  the  value  of  g  required  for 
stability  increases  very  slowly  with  A  .  These  results  are  in  qualitative 
agreement  with  the  results  shown  in  the  Appendix,  which  show  that  this  mode  can 
be  stable  for  positive  A  . 
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9 

Fig.  4  -  Complex  growth  rate  for  odd  5  resistive  interchange  (tearing)  mode  as  a 

* 

function  of  g,  with  ■  -  .06,  S  -  0.7,  -  0.167,  D  ■  0.14. 
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Fig.  8  “  Eigenfunctions  for  g 


Complex  growth  rate  for  even  £  resistive  Interchange  mode  (*.  *  0)  as  a 
function  of  g  with  LA  *  0,  S  *  0.7,  r  S  •  0.0512,  D  «  0.14. 


V.  Conclusions 


We  have  analyzed  Che  stabiliCy  of  tearing  and  resistive  interchange  modes  in 

cylindrical  geometry  for  parameters  relevant  to  current  reversed  field  pinch 

experiments.  Our  principal  conclusions  are  that  these  modes  are  strongly 

stabilized  by  the  inclusion  of  the  Hall  terms  in  Ohm's  law,  the  retention  of 

plasma  compression,  and  exact  treatment  of  field  line  curvature.  The  same 

results  are  not  obtained  from  a  slab  model  which  includes  Hall  terms  but  models 

field  line  curvature  with  a  fictitious  gravity. 

In  the  derivation  of  the  governing  equations  we  made  two  approximations 

which  if  relaxed  may  modify  our  results.  The  first  of  these  was  the  neglect  of 

the  equilibrium  electron  temperature  gradient  and  electron  thermal 

conductivity.  Analysis  of  tearing  modes  in  slab  geometry  has  indicated  a  strong 

21  22 

stabilizing  effect  associated  with  the  equilibrium  temperature  gradient.  ’ 

The  second  approximation  was  <<  Te>  which  allowed  us  to  neglect  gyro- 
viscosity,  and  the  gyro-heat  flux.  In  addition  we  often  assumed  <  =  0.  These 

approximations  have  the  effect  of  eliminating  from  considerations  localized  drift 

23 

modes  which  are  driven  by  the  ion  temperature  gradient  and  drift  modes  which 
are  driven  unstable  by  unfavorable  ion  curvature  drifts. 

Nevertheless,  we  have  seen  that  the  modes  which  evolve  from  the  interchange 
and  tearing  modes  of  a  collisional  plasma  become  stable  as  the  collisionality  is 
decreased  and  the  Hall  terms  become  important. 
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APPENDIX 


In  this  appendix  we  examine  the  approximate  analytic  solution  of  Eqs.  (16)  - 

(18)  for  modes  that  couple  to  the  outer  region  for  the  case  in  which 
3 

g  6  »  1.  The  modes  that  we  consider  have  even  ¥  and  odd  £  and  T.  We  recall  that 

modes  with  the  opposite  parity  were  treated  in  Sec.  Ill,  where  it  was  found  that 

3  2 

the  constant  ¥  approximation  is  valid  provided  g  8  «  1.  We  will  assume  again 

that  the  constant  ¥  approximation  is  valid,  and  thus,  our  results  are  restricted 
to  low  values  of  6. 

In  the  inner  region  ¥  consists  of  a  constant  f  plus  a  correction  ¥^.  We 
integrate  the  second  derivative  of  the  correction  given  by  Eq.  (16)  over  the 
inner  region  and  match  to  the  outer  solution  ¥  —  f  (1  +  1/2  (XlL^A  )  obtaining. 


r  dX  dV/dX2  1  =  r  4^  (Q2  d2£/dX2  +  T).  (Al) 

— i  r  o  -»  \ 


2  2 

Eliminating  d  ¥/dX  in  Eqs.  (16)  and  (17),  replacing  ¥  by  its  constant  value,  and 
dropping  terms  which  are  small  in  the  limit  D  -  8  «  1,  Eqs.  (16)  -  (18)  become 


Q* 


K  +  1(3*>  x¥  +  Qx2e  +  (*ir  x2  ■  0*] 


(A2) 


and 


$  ■  +  »(s  -  ft* 

UA  y  C  C 


21Q*  2  A 

+  Yj-  *<T  “  °e>  -  21Q  «  • 

c  dX 


X¥ 


(A3) 


The  mode  growth  rate  Q  is  determined  by  solving  the  inhomogeneous  equations 
(A2)  and  (A3)  for  T  and  £,  and  inserting  the  solutions  in  Eq.  (Al). 


3 

To  solve  Eqs.  (A2)  and  (A3)  we  assume  g  0  »  1  and  recall  the  results  of 

3 

Sec.  III.  It  was  found  that  in  the  limit  g  S  »  1  two  types  of  localized  mode 
were  found  for  the  homogeneous  system.  The  first  type  of  mode  was  the  drift  wave 
which  had  the  scalings 


X  -  gS1/2,  Q  ~  Q*  -  gB,  DC  -  T. 


Assuming  these  scalings  in  Eqs.  (A2)  and  (A3)  the  dominant  terms  are 


2  ^*2 

(Q  +  iQ*)  XY  +  QX  C  +  i  r  X  T  «  0 
*  o  u 


(A4) 


and 


Q(4  +  r6>  t  +  <!(s  -t2*)  « 

Q  c  c 

2iQ*  7  7  7 

+  TT  -  2iQ  8  d^C/dX^  =  0.  (A3) 


-1  -  1/2 

From  Eqs.  (A4)  and  (A5)  we  see  C  ~  8  B  Thus,  we  can  estimate  the 

contribution  from  the  drift  wave  solution  to  the  integral  in  Eq.  (Al)  to  be  of 

.  „l/2  ~  1 

order  S  g 


The  second  type  of  mode  that  was  found  in  Sec.  Ill  was  a  highly  localized 

-  1/2  3  2 

density  perturbation  which  had  the  scaling  X  «»  g  *  Q  ~  Q*  ~  gB,  T  ~  g  B  C* 
The  dominant  terms  in  Eqs.  (A2)  and  (A3)  for  this  mode  are 


K  [(Q  +  iQ*)  X*  +  (~  x2  -  i)t] , 


(A6) 


and 
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(A7) 


$  -  rs  <« +  ‘“V  t  - 21Q2g  4  • 

dX  c  dX 

1/2 

From  Eqs.  (A6)  and  (A7)  we  can  estimate  T  ~  g  g^.  Consequently,  the 

1/2 

contribution  to  the  integral  in  Eq.  (Al)  will  be  of  order  g  gY^. 

On  the  basis  of  the  previous  estimates  we  conclude  that  in  the 
3 

limit  g  g  »  1  the  integral  in  Eq.  (Al)  is  dominated  by  the  contributions  coming 

from  those  values  of  X  corresponding  to  the  highly  localized  mode,  i.e. 

-  i/2  '  1/2 

X  -  g  .If  one  assumes  -  g  g  then  the  dispersion  relation  is  obtained 

by  solving  Eqs.  (A6)  and  (A7)  for  T  and  C»  and  inserting  the  solutions  in  the 
integral  in  Eq.  (Al). 

Defining  the  normalized  variables  4>,Y,  and  a^,  where 


i(Q  +  1Q*)S3/S 


4  8 


1/2 


♦  . 


1  1/4  -  1/ 2y 

X  =  ^  ®  8  *> 


1/2  1Q*k 

ai  =~2T  ^r\  +  T~T+  ig)' 

c  c 


Equations  (A6)  and  (A7)  can  be  cast  in  the  form  of  an  inhomogeneous  Weber 
equation, 


(A8) 


The  Integral  in  Eq.  (Al)  then  becomes 


-  j  SlMg-  1/2(Q  +  IQ*)  O Y(1  -  \  Y*(Y)).  (A9) 
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Equation  (A8)  may  be  solved  using  the  Fourier  transform  methods  discussed  in  Sec 
III.  The  resulting  dispersion  relation  obtained  form  Eq.  (A9)  is 


21/2*  S1/A 
81/2 


(Q  +  iQ*) 


(A10) 


Equation  (A10)  possesses  an  infinite  number  of  solutions  for  the  growth  rate 

Q  which  are  associated  with  the  poles  of  the  gamma  functions.  These  solutions 

are  the  odd  counterparts  of  the  localized  solutions  discussed  in  Sec.  Ill  and 

have  a  dispersion  relation  similar  to  Eq.  (50).  An  additional  solution  of  Eq. 

-  ^ 

(A10)  is  present  due  to  the  factor  (Q  +  iQ*).  If  A  L  is  small  we  find  Q  — 

*  r 

-  iQ^  +  6Q  where 


1/2 


6Q  =  LrA  |T(^  +  lvj)l  2(coshirv ^  +  i  sinhnv^) 

s 


-1 


(All) 


and 


S1/2  f  2D  ,  x 

vi  *  t1  +  r"es  (,c  ~  x))- 

c 

Equation  (All)  indicates  that  this  mode  will  be  unstable  if  A  >0,  with  marginal 

* 

stability  achieved  at  A  =0. 

The  proceeding  analysis  suggests  that  the  tearing  mode  will  be  unstable  for 

-  '  I/O 

any  A  >0.  However,  in  deriving  Eq.  (A10)  we  have  assumed  L_A  ~  g  0  and 

disregarded  the  contribution  to  the  integral  in  Eq.  (Al)  from  the  drift  wave 

1/2  1/2  -1 
region  X  ~  g0  .  This  contribution  was  estimated  to  be  of  order  8  g  and 

-  \/2  -l 

thus,  Eq.  (All)  breaks  down  for  L^A  -  8  g  .  In  this  case  we  assume 
~  1/2  -1/2 

Q  —  iQ^  +  5Q  where  6Q  -  0  g  «  and  the  contribution  from  the  highly 

1/2  -1 

localized  region  scales  as  0  g  .  To  calculate  the  contribution  from  the  drift 
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wave  region  we  insert  Q  —  IQ^  in  Eqs.  (A4)  and  (A5)  and  solve  for  C  and  T. 
Defining  the  normalized  variables  Y,  and  a^,  where 


X  =  (Qg) 


1/2S-  l/4Yf 


5 


iD  g1/2 
2  Q3/2  S1"2 


and 


a  =  -iS1/2 
a2  2  * 


Equations  (A4)  and  (A5)  can  be  rewritten , 


dY 


-  Y. 


The  contribution  to  the  integral  in  Eq.  (Al)  becomes 


L.A  «  -  j  D(QgS 1/2 )”  1/2  /  JY(1  -  i  Y<KY)), 


which  may  be  evaluated  in  a  manner  analogous  to  Eq.  (A9) 


21/2»D  F(4  +  2  a2) 

r  0l/21/20l/4  1 


Q  g  s 


r(t  +  ia2) 


(A12) 


Combining  the  contributions  from  the  drift  wave  and  localized  regions  we  obtain 
the  dispersion  relation, 
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To  assess  the  importance  of  the  drift  wave  region  we  consider  for  simplicity  the 
case  k  =  1,  S  «  1  and  consequently  a^  —  ~  0*  The  expression  for  the  growth 

rate  6Q  reduces  to 


I  D|  A/  3i  D  ,  .  , 

6Q =  “772  exp  i*r  TdT^  +  V 


g  '  reft 
21/21T  S1/4r(y) 


(A14) 


Equation  (A14)  indicates  that  a  threshold  value  of  A  exists  below  which  the  mode 

is  stable,  and  demonstrates  the  stabilizing  effect  on  the  tearing  mode  of 

20  2 1 

coupling  to  the  drift  wave.  *  If  we  examine  the  definitions  of  g  and  we 
see  that  the  threshold  is  independent  of  n  as  n  +  0. 
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